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Abstract 

We analyze the various approaches to construct exchange-correla- 
tion functionals which are able to describe states of definite spin 
multiplicity in the DFT realm and outline the characteristics of pos- 
sible functionals consistent with the Kohn-Sham theory. To achieve 
this goal the unitary group technique is applied to label many-electron 
states of definite total spin and to calculate the corresponding analogs 
of the Roothaan coupling coefficients. The possibility of using range 
separated Coulomb potential of electron-electron interaction for con- 
structing functionals discriminating multiplet states in the d-shells is 
explored and a tentative system of state-specific functionals, covering 
nontrivial correlations in d-shells of transition metal ions, is proposed 
for the Fe^"*" ions. 

1 Introduction 

Although the Density Functional Theory (DFT) based methods of model- 
ing electronic structure of molecules and solids widely proliferate during last 
decades [1-3], the problem of consistent description of transition metal and 
rare earth compounds with open d- and /-shells, respectively, remains a still 
unresolved, challenging problem in this framework [4]. One of the main rea- 
sons for this failure of the DFT is that the multiplet spin/orbital momentum 
states are generally not easily described within the DFT paradigm. The 
source of that intimate "unfriendlieness" of the DFT to the multiplet states 
lays in the " oversymmetry" of the fundamental quantity pertaining in the 
realm of DFT: the one-electron density. As it has been demonstrated many 
times, states of different total spin and/or spatial symmetry may produce 
equal one-electron densities. The complication arising from this is the im- 
possibility to distinguish the nature of the ground state on the basis of the 
total density only: although only one of say two functions represents the 
ground state, i.e. the exact energies of the involved states may be differ- 
ent, they turn out to be degenerate in the DFT context. In other terms, if 
the same densities are fed to the "universal" density functional implied by 
the DFT, it is going to produce the same value of the electronic energy for 
states whose exact energies are different. Of course the latter remark may 
be opposed by noting that the "universal" functional is going to output the 
ground state energy only, but in this case it is not clear how other important 
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pieces of information concerning the nature of this ground state (e.g. its spin 
multiphcity) can be extracted from such an answer. 

This situation certainly requires some clarification which is addressed in 
the present paper. In order to do so we give below a brief description of 
relevant elements of the electronic structure theory (Section [2]). Then we 
consider an archetypical example of the problems encountered by the DFT 
while trying to reproduce correct spin properties of many-electron systems 
(Section [3]). Then we propose a general scheme allowing to include states of 
definite spin in the DFT theory (Section H]). This, however, does not solve 
the problem of the multiple states in the open d- and /-shells of the transi- 
tion and rare earth ions. For this end we explore in Section E] the possibility 
to circumvent these problems with use of the short/long range separation of 
Coulomb interaction between electrons and propose in Section [5] some con- 
ceivable state-dependent definition of exchange-correlation functionals capa- 
ble to reproduce the energies of nontrivially correlated many-electronic states 
in the ci-shell of the Fe^"*" ions. 

2 Theoretical background 
2.1 Electronic distribution. 

The main idea of the DFT is to reduce the description of entire electronic 
structure to a single quantity: the one-electron density - the diagonal part of 
the one-electron density matrix. The possibility of such a reduction is proven 
by the Hohenberg-Kohn theorems [5] which state an existence of a universal 
one-to-one correspondence between one-electron external potential and the 
one-electron density in that sense that not only the one-electron potential 
acting upon a given number of electrons uniquely defines the ground state of 
such a system i.e. its wave function and thus the one-electron density, but 
also that for each given density integrating to a given number of electrons N 
a one-electron potential yielding that given density can be uniquely defined 
from the density. The "density only" formulation of the electronic structure 
problem, even if it is practically achieved, leaves unanswered an important 
question of the nature of the ground state thus obtained e.g. about its total 
spin (or other symmetry features). 

Incidentally, the symmetry properties of quantum states, like total spin, 
are easier expressed in terms of wave functions (see Ref. [6]) so it would be 
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practical to consider tentative relation between the wave function and den- 
sity only pictures of the electronic structure. The required relation can be 
established with use of the reduced one- and two-electron density matrices 
as much simpler objects than the wave functions, providing equivalent de- 
scription of electronic structure. The reduced density matrices respectively 
depend on two {x,x') and four ( o) coordinates: 



X "^{x' ,X2, ■ . ■ ,Xiy)dX2 ■ ■ ■ dXjy, 
p^'^\xiX2; X[x'2) = J ^^(xi, X2, X3, . . . Xn) X 

X '^{x[,X2,X3, . . . ,XN)dX3 . . .dx^, 



(1) 



where the composite electronic coordinate x represents a pair (r, s) of three 
dimensional radius vector r of an electron and of the discrete variable s 
taking either of the two allowed values ±h/2 of the projection of electronic 
spin. The transition to the description in terms of reduced density matrices 
is itself a significant simplification (although being absolutely exact). The 
one-electron density implied by the DFT theory appears then as a result of 
further reduction of eq. ([1]) : 

pir) = J2p^'\^s;rs) (2) 

s 

Thus according to the DFT paradigm the one-electron density which depends 
on one spatial radius- vector must be able to serve as an equivalent substitute 
to the A^-electronic wave function dependent on N radius vectors and more 
spin projections of all electrons involved. The obvious loss of information 
which takes place by going from eq. ([1]) to eq. ([2]) - we remind that going 
from the wave function X2, Xs, . . . xat) to the reduced density matrices 

by eq. ([1]) does not produce any loss of at least important information - must 
be compensated by the "universal" and "exact" density functional which is 
in general unknown. 



2.2 Electronic energy. 

Leaving aside the "ideal" DFT using the unknown "universal" and "exact" 
functional of the density eq. ([2]) and turning to pragmatic methods pertain- 
ing to the DFT realm needs some approximate expressions for the energy 
presented as a functional of the density eq. ([2]) only. In the wave function 
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and in the equivalent reduced density matrix formulations the energy has the 
form: 

E = {t) + {VUm)) + {Vee)+Vnn{m). (3) 

where {R} stands for the set of radius- vectors of all nuclei inducing the elec- 
trostatic potential external to the electrons of the system. In the coordinate 
representation the above averages acquire familiar appearance: 



t) = -^E j A'pW(rs;r's)t/r 



l^ffp(W£^^,^,^, (4) 



2 ss' - - 

K„({R}) = ^E^^^'*here 



where the expressions eq. (jlj) are assumed to be specific for a given geometry 
{R} and for an electronic state described by the iV-electronic wave function 
^ = \E'(xi,...,XAr) employed to define the density matrices eq. ([T]). The first 
row in eq. (jl]) is the kinetic energy of electrons, the second row is the energy 
of Coulomb attraction of electrons to nuclei, the third row is the energy of 
interelectronic repulsion; the fourth one is the energy of Coulomb repulsion 
of the nuclei, which does not depend on the electronic density/wavefunction. 

In the above expressions eqs ([3]), (jl]) only the average of the nuclear poten- 
tial Vne is exactly a functional of the required form: that of the one-electron 
density eq. ([2]). All other terms in eqs ([3]), (jl]) require further consideration. 
It applies similarly to the remaining one- and two-electron contributions to 
the energy. As for the one-electron term, the kinetic energy requires knowl- 
edge of the one-electron density matrix eq. ([1]) rather than its diagonal part 
eq. although effectively in a narrow range of spatial separations r — r' 
which must be sufficient to determine the second derivative. The attempts 
to avoid this bottleneck and to obtain pragmatic DFT methods brought 
Kohn and Sham [7] to their famous orbital construct which allowed them to 
express the kinetic energy in terms of some single-determinant wave function 
yielding by definition the required (exact) one-electron density. Then the 
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kinetic energy is calculated as one of the system of non-interacting electrons 
described by a single determinant built of KS orbitals. 



2.3 Electronic density and electronic energy decompo- 
sitions. 

While the Kohn-Sham construct offers an efficient technique to handle the 
difficult kinetic energy problem and provide a very good first approximation 
to it, the representation of the electron-electron interaction energy in terms 
of the one-electron density (and possibly further parameters derived from 
the KS determinant) remains the central problem on modern density func- 
tional theory. Generally, calculating the Coulomb electron-electron energy 
(3-rd row of eq. (jl])) requires knowledge of the two-electron density matrix. 
According to [8] it decomposes: 



p(^)(Xi,X2;Xi,X2) = ^ 



p(^)(xi;xi) p(^)(x2;a;'i' 



- x(a;i,X2;x'i,X2), (5) 



where the first (determinant) term corresponds to the part of the two-particle 
density matrix which can be accounted for even in the independent electrons 
approximation. The second term in eq. ([5]) - the cumulant of the two-particle 
density matrix - is responsible for deviation of electrons' behavior from the in- 
dependent electron model, i.e. for their Coulomb correlations. The Coulomb 
interaction of electrons eq. (jl]) can be decomposed to contributions associ- 
ated to the terms of the above two-particle density matrix decomposition eq. 
0: 

{Vee) = Eh + E,c] (6) 

Exc = Ex + Ec. 

by singling out first the "classical" part of the Coulomb interaction energy -~ 
the Hartree energy: 

1 ^ ^ ,W(r..r.)p("(rV.r'.-) _ 1 f pWpM ,,,,, p, 
2^ J r-r' 2 7 r-r' ^ ^ 

ss' 
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and then the exchange and correlation energies: 



E. 



ss' 



drdv' 



(8) 



X 



r — r' 



x(rs, r's'; rs, v's') 



(9) 



r — r' 



whose definitions eqs ([8]), ([9]) are given respectively in terms of the off- 
diagonal part of one-electron density matrix p^^)(r's,rs) and of the two- 
electron density matrix cumulant x{^s, r's'; rs, r's') - the difference between 
the exact two-electron density matrix and its Hartree-Fock (self-consistent 
field) estimate. 

While the definition of the Hartree-energy is unique, and constitutes to- 
gether with the nuclear-electron repulsion energy the part of the total energy 
that can be written straightforwardly as a simple analytic functional of the 
one-particle density, the exchange and correlation energies are defined in 
quantum chemistry and in DFT in different ways. As far as the exchange 
energy is concerned, one should remark that the one-particle density matrix, 
appearing in eq. ([8]) is supposed to be exact. This quantity is not avail- 
able even in exact KS theory, where we have at best the one-particle density 
matrix associated to the single determinant constructed from the exact KS 
orbitals. By consequence, the exact exchange energy in DFT is in general 
not equal to E^. 

It must be observed that the usual definition of the correlation energy in 
quantum chemistry, proposed by Lowdin in Ref. [8] differs from that given in 
Eq. ([H]), which follows rather the suggestion due to Kutzelnigg and Mukher- 
jee [9,10]. This latter definition has the conceptual advantage that it uses 
the quantities entering eqs. ([7]) - ([9]) irrespective to any approximate method 
of calculation of the electronic energy. Some authors, (cf. e.g. Refs. [11,12]) 
argue that pragmatic DFT methods can be considered as approximations to 
the two-electron density matrix cumulant. 

This situation is as well slightly more complicated in conventional Kohn- 
Sham theory, where the correlation energy involves also the difference of the 
exact and KS kinetic energies. However, this kinetic energy contribution can 
be assimilated to a potential energy term the virtue of the adiabatic connec- 
tion procedure, which allows one to write the total correlation energy as an 
average electron-electron interaction over the adiabatic connection path. 
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3 Symmetry non-sensitivity of the density- 
only methods 



3.1 Archetypical example of existing problems 

In order to better understand the problems which appear in the DFT realm 
when trying to describe the correct total spin of a many electronic state 
we consider the simplest system of two electrons occupying spatial orbitals 
I a) and \b) (which can be understood as notation for one-dimensional irre- 
ducible representations of a point group) and forming corresponding singlet 
and triplet states and ^B. The relevant wave functions in the coordinate 
representation are given by: 

^ifilxi, xa) = I (a(ri)6(r2) + 6(ri)a(r2)) {a{si)f3{s2) - /3(si)a(s2)) , 
^3b(xi,X2) = i(a(ri)6(r2) -Kri)a(r2))(a(si)/3(s2) +/5(si)a(s2)), ^ ^ 

both having the zero projection of the total spin. Following the definitions 
of the one-electron density matrices eq. ([1]) the states eq. ffTOj) immediately 
yield exactly the same one-electron density matrix: 

p« = ^ ia*is)a{s')+P*is)Pis')) (a*(r)a(r') + 6* (r)6(r0) (11) 

irrespective to the total spin of these states. This result is well known for 
decades and appears even in textbooks [13]. Obviously the density eq. ([2]) 
which is required by the DFT is as well the same for the two spin states. 

The exact two electron density matrices calculated according to their 
definition eq. from the wave functions eq. ( |TOl) are, however, different: 

i (a*(si)r (S2) T r (si)a*(s2)) {a{s[)(3{s',) T m)a{s',)) x 
(a*(ri)6*(r2) ± 6*(ri)a*(r2)) {a{r[)b{r',) ± b{r[)a{r'^)) 

with the upper sign corresponding to S' = and the lower one to 5"= 1. Com- 
paring the above expression with the decomposition eq. ([5]) one easily sees 
that only the cumulant of the two-electron density matrix can be responsible 
for the distinguishing of the two-electron density matrices for the singlet and 
triplet states. 

The " oversymmetry" of the density (and even of the first order density 
matrix) with respect to the total spin exemplified by eq. ffTTl) is not acci- 
dental, but is a consequence of a very general result (see Ref. [14, 15] and 
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references therein). Even a higher symmetry can be proven [16]. In its mod- 
ern formulation (Theorem 1 of Ref. [14]) it reads: 'The electron density of 
an arbitrary A^-electron system, characterized by the A^-electron wave func- 
tion corresponding to the total spin S and constructed on some orthonormal 
orbital set, does not depend upon the total spin S of the state and always 
preserves the same form as it is for a single-determinant wave function'. The 
proof given in Ref. [14] relies not upon the spin properties themselves rather 
on the manifestation of permutation symmetry of the exact wave function in 
terms of the total spin. We address this issue later in Section 14.2.11 

3.2 Methods proposed to treat coinciding densities 
3.2.1 Multiplet sum method 

The first attempt to get around this problem of coinciding densities in the 
DFT context dates back to the work Ref. [17]. The analysis of problems 
performed there is precisely repeated in the above two-electron two-orbital 
model. The prescription Ref. [17] concerning the way out reads as follows: 
to evaluate correctly the energy of the singlet and triplet states and in 
terms of the quantities which can be obtained with use of single determinant 
wave functions. To do so one has to address the single determinant function 
|aa6/?| which is not a pure spin state, but in fact is a linear combination of 
two above spin states: 

\aab(3\ = {\'B, = O) + \'B, = 0» (12) 

Averaging the Hamiltonian over the linear combination eq. f|T2l) of the pure 
spin states immediately yields: 

^-ECB) + ^-ECB) (13) 

The energy of the triplet state entering the above combination can be inde- 
pendently extracted from another single determinant wave function: |aa6a| 
corresponding to the component of the triplet with the spin projection -|-1. 
Thus one can express the energy of the (non-single-determinant) singlet state 
linearly combining the averages of the Hamiltonian over the single determi- 
nant states of which, however, one belongs to the spin projection -|-1. Ob- 
viously the above move was only possible because the off-diagonal matrix 



9 



element of the Hamiltonian between the singlet and triplet contributions to 
the determinant of interest vanishes due to the spin symmetry. The different 
expressions for E(^B) and E{^B) thus obtained are then treated as required 
distinct energy functionals to be used to calculate the energy respectively 
for the singlet and triplet states possessing the same one-electron density. 
It is instructive to check (and in this simple case it can be done by direct 
evaluation) where the difference between the energy expressions comes from. 
Inserting the one-electron density matrix eq. (fTT!) which is the same for both 
spin states in the definitions of the Hartree and exchange energies eqs. ([7]), 
(IH]) yields for the both spin states equal Hartree and exchange contributions: 



One can see that (i) the self interaction terms in the Hartree contribution 
are precisely cancelled by the corresponding terms in the exchange part; (ii) 
at the same time, obviously, there is no other source where the difference 
between the spin state energies could come from except the cumulant of the 
two-electron density matrix and thus the correlation energy as defined by eq. 
(ini) is responsible for the difference in the resulting expressions: 



On the other hand one may notice that the classification of the energy con- 
tributions as exchange or correlation ones by eqs. ([8]), ([9]) is in some way 
arbitrary as well. Indeed, for the above model the energy of the triplet state 
with the spin projection -|-1 is exactly the sum of the Hartree and exchange 
contributions since the latter state is represented by a single determinant 
wave function for which the cumulant precisely vanishes. However, the equal 
energy for the triplet state with the zero spin projection breaks down differ- 
ently: into the Hartree, exchange, and correlation contributions, where the 
Hartree contribution is the same as in the case of the spin projection but 
only the sum of the exchange and correlation contributions is the same for 
different values of Sz- 

The above way leading to the energy expressions for different spin states is 
not completely satisfactory (and not clearly generalizable) : although formally 
the results can be treated as functionals of the density the difference of the 
two energy expressions is obtained by a kind of trick. Referring to the triplet 



Hartree \ [{aa\aa) -\- {aa\bh) + {hh\aa) + {hh\bh)] ; 
exchange ^ [(aa|aa) + {ab\ba) + {ba\ab) + {bb\bb)] . 



(14) 



E(^B) = {aa\bb) + {ab\ba) 
E{^B) = {aa\bb) - {ab\ba) 



(15) 
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component with S'^ = +1 in the derivation of the multiplet energy looks out 
as an ahen element (in fact the energy is uniquely determined by the spatial 
multipliers in the wave functions eq. ffTOl) - without any reference to the spin 
components at all). This strange element of the derivation appeared in order 
to compensate somehow the element of the general theory which is missing 
in the DFT - the cumulant of the two-electron density matrix. Despite 
this criticism the result of the derivation is very transparent: it reduces to 
deriving according to McWeeny's notice in Ref. [18] 'of a particular type of 
energy expression — irrespective of the nature of the wavefunction', namely 
one - linear in the Coulomb and exchange two-electron integrals over the 
involved orbitals. 

Further development of this approach is based on the assumption that it is 
always (or at least for unspecifically wide class of cases) possible to express 
the energy of a pure spin multiplet state allowable for a given number of 
electrons/orbitals as a linear combination: 



where Ei are the diagonal matrix elements of the energy operator taken 
with respect to all necessary Slater determinants (numbered by i). The one- 
electron density (matrices) corresponding to these determinants are different 
and the whole scheme becomes workable provided the set of the coefficients 
(weights) wf^ exists and they are uniquely defined by the spin and sym- 
metry quantum numbers F and S and other quantum numbers n serving to 
distinguish potentially existing states with equal F and S. Apparently only a 
restricted number of examples of such functional forms is known. The reason 
is quite simple and the above consideration allows to single out the range of 
cases where the derivation analogous to that of eq. f|T2l) can be performed. It 
applies if the total spin allows to completely distinguish the electronic states. 
In this case the energy of the state of the highest available spin S'max can be 
expressed through a single determinant function (the cumulant is vanishing) 
with the highest available projection of the total spin. Then this result can 
be used to evaluate the energy of the state with S'max — 1, etc. The recipe 
immediately fails as soon as multiple states of the same total spin appear in 
the system. This is however the everyday life, so in what follows we switch to 
considering further possibilities of constructing the energy functionals useful 
in this situation. 




(16) 
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3.2.2 Restricted open shell KS (ROKS) method 

The situation with reproducing total spin dependence of the energy as it 
appears in the DFT context is by no means unique: the same problem arises 
in the Hartree-Fock-Roothaan (HFR) context since the latter lacks any ad- 
equate representation of the cumulant of the two-electron density matrix as 
well. Within the "extended" HFR context some ways out have been pro- 
posed. Incidentally, the method of Ref. [17] is precisely the Slater multiplet 
sum method Ref. [19] which migrated from the HFR to the DFT context. 
Another option is the ROHF (restricted open shell Hartree - Fock) method 
whose respective migration resulted in a range of the ROKS (restricted open 
shell Kohn-Sham procedures [20,21]) being the DFT counterpart of the for- 
mer. Despite different appearance they have many common features (and 
we do not address here the methods based on the statistical - ensemble - 
averaging) . 

The derivation of the ROHF (or equivalently 'old MC SCF ' - see below) 
bases on the general expression of the form: 



where C^^^ are the expansion coefficients of the eigenfunction '^nrs{xi,X2, ^3, . . . xjy) 
of the many-electron hamiltonian over some appropriate basis states ^i{xi,X2, x^, . . . 
[e.g. the Slater determinants, see, however, below). In this case no alien 
states of wrong spin projection may appear. On the other hand the contri- 
bution of the off-diagonal elements Hij to the energy may be nontrivial (in 
contrast with eq. f|T6l) ). The knowledge of the expansion coefficients C'^^^ in 
general requires diagonalization of the Hamiltonian matrix making the ex- 
pansion coefficients and thus the energy itself some sophisticated irrational 
function of the Hamiltonian matrix elements including two-electron integrals. 
It was Roothaan [22] who first noticed that certain states "^nvs of atoms and 
linear molecules, even those requiring many-determinant (multi-reference, 
multi-configurational) wave functions, yield energy expressions which are 
linear with respect to two-electron integrals {ii\jj) and (respectively 
Coulomb and exchange ones). It is only possible if the wave function ex- 
pansion coefficients C^^^ in eq. f|T7j) can be determined on the symmetry 
grounds i.e. without nontrivial diagonalization. In this case there is no 
need that the off-diagonal elements Hij which are linear expressions in the 
two-electron integrals and thus give a linear contribution to the energy func- 




(17) 
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tional disappear as required by eq. (fTBI) . Only the possibility to have the 
expansion coefficients independent on the specific values of the Hamil- 
tonian matrix elements is the true prerequisite for obtaining the expressions 
for the energy of the required (linear) form. Nevertheless the number of 
cases when the described procedure was possible in fact reduces to the 
states of atoms and vr" and 5" states of linear molecules. Similarly the ROKS 
sheme proposed in Ref. [21] and representing a migration of the Roothaan's 
reasoning to the DFT context allowed to obtain the functional forms for the 
same set of states: p", vr"-, and (5"". Thus the forecast of the year 1960 due 
to Roothaan [22]: 'It is a relatively simple matter to extend the open-shell 
theory just presented in such a way that other important classes of atomic 
states can be accommodated, as for instance, the configurations for the 
transition elements. We postpone such generalizations for the present, and 
include whatever new treatments may be necessary with the actual applica- 
tions planned for the future' never became true and the states generally 
cannot be squeezed in the ROHF/ROKS scheme. 

Under other angle of view, validity of the Roothaan or similar schemes 
means that the cumulant of the two-electron density matrix can be in some 
particular case recovered by symmetry based manipulations. In the cases 
considered by Roothaan himself and recently used in the DFT context in 
Ref. [21] the possibility of obtaining closed expressions for the energy func- 
tional in terms of two-electron integrals over orbitals involved is stipulated 
by additional symmetry of the system (in the chemical context it goes about 
additional symmetry group G with irreducible representations F, where G = 
5*0(3) for an atom G = SO {2) for a linear molecule, and may be some 
point group for other molecules) which allows to figure out the expansion 
coefficients C^^^ . It is clear that for an overwhelming majority of cases it 
is impossible to find any nontrivial symmetry group G ^ C\ which prede- 
fines a restricted character of any Roothaan-like treatment. It is thus our 
next purpose is to expore other possibilities of designing energy functional 
distinguishing the states of different spin multiplicities in the DFT context. 
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4 Spin and unitary symmetry of the electronic 
wave function 



As wc mentioned previously any reference to the spin projections through- 
out the derivation of the energy expressions for the two-electron two-orbital 
model looks out as an alien element. The ultimate reason for that is that 
the non-relativistic Hamiltonian does not depend on spin variables at all and 
the energy itself as well as the differences in its form for different spin states 
originates solely from the spatial multiplier of the many-electronic wave func- 
tion (spatial function). The idea to restrict the entire consideration by those 
spatial functions persists almost from the beginning of the quantum chem- 
istry and is known as "spin-free quantum chemistry" [23]. It can be given 
different formulations of which we use one based on the use of the unitary 
group (see Ref. [18]). We briefly remind its basic facts in the Appendix. 

4.1 Unitary symmetry of the spatial multiplier. 

The construct using the permutational symmetry of the spatial part of the 
wave function had been used for developing the so called generalized Hartree- 
Fock procedure [24] which had numerous descendants (see e.g. [25,26]). They 
basically performed the task of presenting the energy in the HFR-like form: 
linear with respect to Coulomb and exchange integrals over the orbitals in- 
volved with the coefficients dependent on the permutational symmetry of the 
spatial part of the wave function and thus on the total spin. The permuta- 
tional symmetry, however, addresses the many-electron wave functions in the 
coordinate representation which is of restricted use in quantum chemistry. 
By contrast the wave functions actually used are those in the representation 
of the occupation numbers of the orbitals involved. For that reason it is more 
practical to switch to labeling of the many-electron functions by irreducible 
representations of the unitary group which are closely related to those of the 
Sjsi group. The corresponding construct is described in the Appendix. 

4.2 Physical quantities in terms of unitary group. 

Going to the representation of the U{M) group has that advantage that 
it allows to easily write down the energy of many electron states. This is 
done as follows: for each Young pattern T one can construct the set of 
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generators {ij = 1 -^ M) of the group U (M) acting in the space of the 
irreducible representation T = T{M,N,S) whose matrix elements between 
the tableaux v and v' can be calculated irrespective to the physical nature 
of the system described. The set of generators completely defines the action 
of the group U{M) in the irreducible subspace of its tensors of the rank N 
with the permutational/spin symmetry stipulated by the Young pattern T. 

The diagonal generators are diagonal in the basis of Young tableaux 
and their matrix elements are equal to the occupation number (rij = 2, 1, 0) 
of the i-th orbital in the Young tableau Tv : 

{Tv |E,^| Tv') = 5^^, {^l^tv = ^^^'^i 

By contrast off-diagonal generators E^ (raising ones if i > j and lower- 
ing ones if j > i) have non- vanishing matrix elements (Tv |E^| Tv'^ if the 
tableau v' contains at least one orbital symbol j whereas the tableau v con- 
tains one less orbital symbol j than v' and one more orbital symbol i than 
it. From this it follows that the off-diagonal generators E^(i ^ j) have no 
non- vanishing diagonal matrix elements. 

4.2.1 One-electron density in the unitary group formalism 

The generators E^ are by definition the components of the spatial one- 
electron density operator restricted to the subspace of the states belonging to 
the T pattern (those having transformation properties of the corresponding 
rank tensors with the permutational symmetry stipulated by the Young 
pattern T or equivalently having the total spin prescribed by this pattern): 

a 

where stands for the operator projecting A^-electron wave function to 
the subspace of the functions with the spatial part having the permutational 
symmetry of the irreducible representation T of either or U{M) groups. 
It is remarkable to note that the Young tableaux states have an important 
property similar to that of the Slater determinants: the one-electron density 
matrices generated from such states are diagonal in the basis of the orbitals 
involved in their construction. 

With use of this construct one can easily check the validity of the Kaplan's 
Theorem 1. Indeed, for whatever Young tableau Tv the one electron density 
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pertinent to the corresponding A^-electron state reads: 
Pr,{r,r') = ^ p« (rs, r's) = (Tz; |V^+(rs)V^(r'5) | T^;) (18) 

s s 
(7 s ij 

ij 
i 

which in turn does not depend on the permutation symmetry labels Tv, 
which is the only connection to the total spin. Thus even the spatial density 
matrix (not only the density) is permutation/spin independent as stated in 
Ref. [14]. 

4.2.2 Energy in the unitciry group formalism 

Further development is based on the possibihty to express the blocks of 
the Hamiltonian matrix pertaining to N electrons in M orbitals with total 
spin S through the generators E^, with T = T{M, N, S). The required 
representation reads Ref. [18]: 

H = ©H^ (19) 

T 

= J2h,,El + ^Y.(tm{ElEl-Sj,B'^) (20) 

ij ijkl 

The matrix elements hij are the sums of the respective matrix elements of 
the kinetic energy Te of electrons and of the external Coulomb potential Vne', 
the quantities {ij\kl) - the two-electron matrix elements of the Coulomb 

interactions. 

For each of the states \Tv) represented by the Young tableau with the 
Young pattern T and the fiUing v (this information suffice to define the 
spatial part of the A'"-electron wave function) the expectation value of the 
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energy reads: 



EiTv) = h., (Ej)^^ + I {{Ej^El - S,,El))^^ (21) 



2 

ij ijkl 

For the one-electron contribution to the energy one gets: 

E (Ej>T. = E (^)r. = E ^^^^^ 

ij i i 

and the Coulomb interaction of electrons is expressed through the Coulomb 
and exchange integrals with respect to the orbitals involved in the construc- 
tion of the states represented by the Young tableaux: 



Hartree i E (E,^EJ^>^„ + 

i^j i 



The Young tableau states Tv are the eigenstates of the diagonal generators 
E,^. For that reason the Hartree contribution to the energy can be rewritten: 

Hartree i E {^I)r. (^Dtv = \ E i^^Jj) n^n^ (23) 

^j ij 

in terms of the products of the one-electron densities. From this we see 
that the Hartree part of the Coulomb energy is uniquely defined by the 
occupation numbers of the spatial orbitals i.e. only by the spatial density 
in the representation of orbitals. We see that as in the other representations 
the Hartree term is contaminated by the self-interaction of electrons and 
that the principal effect for which the true exchange term is responsible in 
the HFR context - the avoiding of the self interaction - is guaranteed by the 
specific form of the coefficient at the integrals of the type which must 

be absorbed by the exchange contribution to the energy. The averages of the 
off-diagonal generators' products entering the expression eq. (12^ are not 
however uniquely defined either by the occupation numbers of the orbitals 
in the tableau Tv or by the total spin, prescribed by the pattern T. They 
depend also on the mutual positions of the orbital symbols in the tableau. 
This is precisely the result obtained many years ago in Ref. [27] under the 
name of the spin-free self consistent field theory. 
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/^Frorn the ROHF (old MCSCF) point of view the result eqs. ([21]), 
can be considered as a recipe of obtaining the coupling coefficients and 
bij at the Coulomb and exchange integrals in the ROHF expressions for the 
energy which incidentally acquire the Tv dependence: 

aj/^ = {Ej^El-6,,El)^^ (24) 

= {EM-^r. (25) 

Turning back to expressions eq . (fT7|) one can say that constructing the 
spatial Young tableaux states provide the expansion coefficients C^" for the 
respective linear combinations of the A^-electron Slater determinants, yielding 
the total spin specified by the Young pattern. These coefficients are derived 
by purely symmetry reasons and do not depend on the matrix elements of 
the Hamiltonian thus satisfying the requirement of " universality" . On the 
other hand it is obvious that specifying the total spin only does not suffice to 
specify the electronic state. The procedure implied by eq. fl2T!) provides for 
each allowable set of M, N, S the whole bunch of energy expressions labeled 
by the rows v of the irreducible representation T = T(M, N, S). 



4.3 Multiplet sum method from the unitary perspec- 
tive 

The first usage of the above formalism is to repeat the success of the MSM in 
case of two electrons in two orbitals without addressing explicitly the foreign 
component of the triplet state with = +1. Indeed, the spatial parts of 
the multiplet states in eq. (fTO|l are equivalently represented as the Young 
tableaux states: 



'B 



Two electrons in two orbitals form only one spatial function for the spin 
triplet state, but in addition to one given above two more functions compat- 
ible with the spin singlet state: 

> 



a 


a 


>> 


b 


b 



are available. Three spatial functions compatible with the spin singlet state 
together form a basis of the three-dimensional irreducible representation of 
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the group U{2) corresponding to the total spin 0. The single spatial function 
compatible with the spin triplet state spans the one- dimensional irreducible 
representation of the group U{2). The Young pattern label T here can be 
replaced by indicating the total spin only. Then the generator E^^^ = 0, but 
for S = one has: 



a 


a 




a 


b 


) = V2 = { 


h 


h 


^ab 


a 


b 



These values suffice to perform the matrix multiplication of the generators 
Ef^°E£^° in the general expressions eq. ( l22l) . so that we obtain for the 
contribution of the average interaction to the energy: 

E(^B) = {aa\bb) + {ab\ba) 
e\^B) = laa\bb) - lab\ba) 

as it should be. We see that the archetypical result is reproduced within 
the Young tableaux technique without addressing the component of the spin 
multiplet with a foreign value of the spin projection. Also the self-interaction 
contamination is removed automatically. 

4.4 DFT implications 

All above treatment was not in any way related to the DFT realm. The 
possibility of establishing such a relation can be based on the recognition of 
the fact that the symmetry (in particular the total spin) dependence must 
be extraneously introduced into DFT considerations [28,29] analogously to 
the treatment by Filatov and Shaik Ref. [21]. The unitary group formalism 
allows us to conclude that for a given set of consistent values of M, N, S one 
arrives to the family of functionals labelled by the rows v of the irreducible 
representation T = T{M,N,S) of the U{M) group. The spin symmetry 
features of these functionals are condensed in the (EJjEJ^ — E^)^^ (or ajj^, 
bjj") coefficients given above. 

The energy matrix elements reflecting specific features of the system can 
be easily figured out. The coefficients ajj" for the Coulomb integrals 
which define the Hartree part of the Coulomb energy are known, but they 
are of no practical use in the DFT context, where the Hartree part of the 
interaction is determined directly from the electron density. Relatively prob- 
lematic (in the DFT context) is to decide where the energy matrix element 
to be combined with the coupling coefficients (E^E^ — E,^)^^ (exchange) 
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and — (E^)^^ (self-interaction) has to come from. This choice must be com- 
patible with various theoretical settings. First of all we notice that if a 
hybrid functional is used which contains some fraction of the Hartree-Fock 
exchange the latter must be modified accordingly so that the corresponding 
integrals over the Kohn-Sham orbitals be included with the correct 
coefficients (E^EJ^ — E^)^^. The same apphes to the integrals which 
together with coefficients — (E^)^^ will take care about some fraction of 
self-interaction. 

Further concerns are related with the treatment of the nontrivial parts 
of the exchange-correlation functionals within the Tv numbering of the spin 
(permutation) states. This can be solved on the basis of certain consistency 
requirements. Among possible consistency requirements the most natural 
is that with the TDDFT. The TDDFT approximation is equivalent to con- 
structing the corresponding time evolution of the many-electronic state in 
the basis of single electron excitations (particle-hole pairs) above the KS sin- 
gle determinant wave function. Leaving aside the question of the area of 
applicability of such an approach we notice that it requires an estimate of 
the two-electron integrals coupling between different singly excited determi- 
nants. The interaction appears as second functional derivative of the energy 
with respect to density (first functional derivative of the exchange-correlation 
functional) . In the orbital representation these derivatives acquire the neces- 
sary form of two-electron integrals with the kernels determined by the form 
of the used exchange-correlation functional. On the other hand the 
integrals appear in ROHF and in unitary group formalism for Tv states as 
a consequence of configuration interaction between different Slater determi- 
nants implicitly entering in the expansion of the Young tableau state Tv. 
Thus in order to ensure the compatibility of the corresponding components 
of the theory: the couplings between the elementary excitations in TD-DFT 
and between Slater determinants in expansions of Tv states must be the 
same. Thus they can be expressed through the integral kernels of the inter- 
action (e.g. according to [30]): 

=11 ^:(r)(^,(r)/-(r,rO(^*(rO(^,(r')c^rdr'; 



/-(r,r') 



where (26) 
5p{r') 

Namely these quantities must be inserted in the expressions for the exchange- 
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correlation energies to get these later consistent with the total spin/permutation 
symmetry of the underlying many-electronic ground state. This is also in 
agreement with the way of constructing the coupling operators by Filatov 
and Shaik in their version of ROKS Ref. [21] and a similar procedure can be 
easily designed for the Tv labelled states. 



4.5 Further examples 



As we mentioned many times the spin in general does not suffice to dis- 
tinguish many electronic states with the same one-electron density which 
produces problems in describing corresponding states in the DFT. The only 
example of the usage of the unitary group formalism given so far was how- 
ever the simplest case when the total spin labeling was sufficient. Below we 
briefly exemplify the features one should expect in general case when there 
exist Young tableaux differing by the positions of the orbital symbols in these 
tableaux. In this case one can say that for given M, A^, S uniquely defining 
the irreducible representation T of the group U{M) and for the row v of the 
latter defined by a specific location of the orbital symbols in the tableau a 
" Hartree-Fock-like" energy functional can be written whose electron-electron 
interaction part is given by eq. fl22|) . It can be optimized with respect to 
the expansion coefficients of the involved orbitals over the AO's basis yield- 
ing an effective Fockian matrix whose eigenvectors are precisely the orbitals 
involved in the construction of the Young tableau state in the same manner 
as it is in the ROHF/ROKS. 

It is easy to check that the positions of the orbital indices in the tableaux 
really matter. For example, for two Young tableaux states: 



\Tv) 



a 


b 


c 


d 



|Tt;') 



a 


c 


b 


d 



(27) 



both representing singlet states of four electrons in four orbitals with equal 
one-electron density matrices, the contributions to the energy functionals of 
the form eq fl22l) . proportional to the exchange integrals, respectively, are [31]: 

|Tt;) : {ab\ba) + {cd\dc) — \ [{ac\ca) + {ad\da) + {hc\ch) + {hd\dh)\ 
[Iv') : -lab\ba) - {cd\dc) + \ \{ac\ca) + \ad\da) + \hc\do) + \hd\dh)\ 

(28) 

where the Hartree and the self-interaction correcting terms are omitted for 
brevity. 
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Remarkably enough neither of the expressions eq. 0281) (combined with 
other necessary temrs) yields a lower energy a priori: which one is lower 
depends on the relations between the molecular integrals involved. At this 
point one can return to the qualitative interpretation of the Young tableaux 
with different positions of the orbital symbols as of reflecting different "pair- 
ing schemes" . Indeed, the states in eq. flTTI) can be respectively treated (and 
this is in accord with the energy expressions eq. fl25]) ) as pairwisely coupling 
electrons in the states a and h and c and d to the singlets and triplets then 
coupling these intermediate states to the final singlet states. 

On the other hand one can easily conclude that for the above pair of 
Young tableaux Tv and Tv' for which = rij = 1 and the difference between 
them is only the positions of the orbital symbols h and c in the tableaux 
eq. fl27|) the operators E^E^ entering as multipliers of the {ii\ii) exchange 
integrals in the exact Hamiltonian yield also an off-diagonal matrix element 
of the Hamiltonian [31]: 

{Tv |H^| Tv') = {{ac\ca) - {ad\da) - {bc\cb) + {bd\db)) 0, (29) 

which shows that in this case a 2x2 diagonalization is required for obtaining 
the electronic energy exactly. Thus the energy functional becomes a square 
root irrational function of the two-electron integrals rather than a linear one. 
Additional symmetry relations may produce the energy expression linear in 
the Coulomb and the exchange integrals (in the above example it suffice that 
exchange integrals {ac\ca)^ {ad\da) and (6c|c6), (6(i|c/6) are pair-wisely equal 
(which makes eq. (129|) be zero) or alternatively that the exchange integrals 
satisfy the equality: 

{ab\ba) -|- [cd\dc) = - [{ac\ca) + {ad\da) + {bc\cb) + {bd\db)] (30) 

which makes the diagonal matrix elements of the Hamiltonian be equal for the 
states represented by eq. (1271) . Both symmetries yield specific forms of the 
2x2 configuration interaction matrix and by this allow the diagonalization 
to be feasible on the purely symmetry grounds). 
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5 State-specific exchange-correlation function- 
als for atomic d-shells 



The above notion of irrationality shows that even the unitary group formal- 
ism does not solve the problem of constructing density functionals for the 
specific correlated states. Although the unitary group formalism allows to 
significantly contract the expansions of the states of the definite total spin (in 
fact the Tv labeled states become single-configuration, albeit each of them 
is a combination of many Slater determinants) the nonlinearity of the energy 
expression with respect to the two-electron integrals hinders constructing the 
symmetry adapted functionals along the lines suggested above. This prob- 
lem manifests itself in the description of the many-electronic states in the 
(i-shells of transition metal ions. Using the unitary goup formalism also in 
this case does not allow to go further than the Roothaan old MC SCF the- 
ory as described in Ref. [22]. There the spin/angular momentum dependent 
coupling coefficients a^j and bij had been introduced ultimately to express 
the cumulant of the two-electron density matrix using symmetry considera- 
tions. They are valid only if the multiplet states can be uniquely obtained 
by applying operators projecting the Young tableau states to the specific 
rows of the irreducible representations of the SO (3) or SO (2) groups (atoms 
and linear molecules, respectively). In these two cases moderately simple 
expressions for the classifying operators (respectively, and L^) in terms 
of the generators E^- can be written and used for constructing the required 
symmetry adapted combinations of the Young tableau states. In the case 
of the atomic p-shells (and molecular tt- and 5-shells) the number of the 
5*0(3) (50(2)) symmetry labels (different values of the orbital momentum 
L) produced by the projection of the Young tableau states to the definite 
states suffice to distinguish all different energies in these shells. In the case 
of atomic p^-states the symmetry 50(3) reduces also the number of inde- 
pendent two-electron integrals (including bothe the Coulomb and exchange 
ones) to only two independent Slater-Condon parameters Fk{pp);k = 0,2. 
This allowed the authors of Ref. [32, 33] to develop state-specific functionals 
for the atomic p^-states. It turns out, however, that for the (i-shells it does 
not suffice for a major part of the atomic electronic terms of the transition 
metal ions [34]. Even in free ions where the multiple terms having the same 
spin and orbital momentum do exist and their energies cannot be expressed 
linearly through the two-electron integrals. Despite the high-symmetry sit- 
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uation of a free atom (ion) which reduces all the two-electron integrals to a 
limited number (three) Slater-Condon parameters Fk\k = 0,2,4 in the free 
ions the energies of the multiplets require 2x2-diagonalization and thus their 
analytical expressions contain square roots (for a handy reference see [34]). 
This moment is crucial - it is not possible to get rid out of the irrational- 
ity (square root) in the expression for the energy by linearly combining the 
parameters of the Hamiltonian. 

The situation clearly becomes less favorable in lower symmetries or in 
larger subshells {e.g. partially filled /-shells) where the terms of the same 
spin and symmetry span the subspaces of dimensionalities higher than two. 
For example, in the octahedral environment the LS states of d^- configu- 
ration span up to seven-dimensional subspaces of many-electronic states [35]. 
Clearly, at an arbitrarily low symmetry the problem of linearly expressing the 
exact energy of many-electronic terms through the Coulomb and exchange 
integrals cannot be solved and obviously the energy of any of such multiple 
terms cannot be expressed as a linear combination of Coulomb and exchange 
integrals. In what follows below we restrict ourselves to the case of atomic d- 
shells and the square root irrationalities in the state-specific expressions for 
the energy trying to squeeze the simplest thinkable irrationality reflecting 
nontrivial correlations in a kind of generalized density functional. 

5.1 The example of Fe^"*" ion 

We concentrate on the free Fe^+ (d^) ion which is an important object in the 
studies of biologically active transition metal complexes and following [34] 
provides a rich system of nontrivially correlated multiple states in its d-shell. 
Namely this kind of behavior is known to systematically evade from any 
DFT-based treatment. The energy expressions of the states in a free Fe^+ 
(d^) ion are given in Table 1. They nontrivially depend on two Slater-Condon 
parameters: F2 and F4. The ground state follows the Hund's rule and for the 
Fe^"*" (d^) ion it is the state. According to the data published Ref. [36] 
the states \.S and j^D are not resolved from the spectra. Also the state 
cannot be present in the spectrum of an even-electron system. Thus we 
exclude three uppermost rows of the Table 1 and flnally arrive to the set of 
data Table 2 which can be used for analysis. 

In order to get an impression of what can be (and should be) possibly 
achieved in terms of describing the multiple states of the d-shells we deter- 
mine parameters F2 and F4 from experimental data. This can be done in a 
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number of ways. The semi-empirical approach is to assume that F2 and F4 
are independent parameters. At the first stage we neglect the correlation and 
take into consideration only the average energies of the ^2^~^^^L states. The 
corresponding set of energies is given in Table 3. These energies are linear in 
the parameters F2 and F4. Applying the standard linear least squares pro- 
cedure yields the experimental "non-correlated" estimate of the parameters 
(in cm~^): 

^exp ^ 1411.0, 

FfP = 120.25, (31) 
FfP/FfP = 11.734. 

The quality of this result can be assessed by the value of mean square devi- 
ation which is 686.80 cm~^ which must be compared with the range of the 
energies described by the model being ca. 45000 cm~^. 

Next step consists in estimating the manifestations of correlations in the 
available data set. The most direct way to do that is to consider the square 
root contributions to the energies of the multiple terms of the same spin and 
symmetry. These come from the diagonalization of the symmetry adapted 
CI matrices. Technically the correlations of that sort are responsible for the 
splitting within the pairs of states of the same spin and symmetry which 
do not have any counterpart in the DFT and describe the nontrivial (non- 
dynamical) part of the correlation. We can see from the Table 4 that the 
correlation splitting between the double states is by approximately 10% un- 
derestimated when calculated with use of the non-correlated experimental 
estimates of the F2^^ and F^^^ parameters eq. (1ST]) . The overall picture as 
coming from the non-correlated estimate can be characterized by its mean 
square deviation 1228.0 cm~^. This fit seems to be improvable by perform- 
ing another (nonlinear) one for the entire set of available excitation energy 
expressions and the corresponding experimental values as given in Table 2. 
The result of this new fit is, however, twofold. The resulting values (cm~^) 
of the parameters F^'''' and F^^^ eq. ( |32|) : 

Fa^P = 1468.92, 

FfP = 113.30, (32) 
FfP/FfP = 12.960, 

which can be qualified as "correlated" experimental ones, produce the mean 
square deviation of 842.37 cm~^ which manifests a significant improvement 
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as compared to analogous usage of the "non-correlated" experimental values 
eq. (13 ip . Meanwhile, although the overall picture is improved the descrip- 
tion of the average multiplet energies is deteriorated as compared to the 
"non-correlated" parameters eq. fl3T|) so that the corresponding mean square 
deviation somewhat increases to the value of 859.37 cm~^. 

The above results deserve thorough attention. First of all we notice fol- 
lowing Refs. [33, 37] that the parameters F2 and F4 are by definition some 
functionals of radial density: 

[ [ [min(ri,r2)] ^ 2/ n 2 2 . > 



Dk J J [max(ri, r2)]'^"'""'^ (33) 
= 1;D2 = 49;L>4 = 441 







where R^{r) is the radial density distribution for the involved atomic d-shell. 
However, according to the Theorem 2 of Ref. [14] for whatever spatial multi- 
plet the one-electron density is spherically symmetric. Thus the quantities 
are the functionals of one-electron density which in the said case have only the 
radial dependence r = |r|. For that reason the energies in Table 1 can be also 
treated as functionals of the one-electron density representing the averages of 
the electron-electron interaction energy for each specific many-electron state 
in the d-shell. When supplied by the relevant one-electron contributions (ex- 
pression for the kinetic energy and that for the electron-nuclear attraction) 
they become the state specific energy functionals 

T[R\t)] + VnAR\r)] + !^li!|^A[i?2(r)] + XC^Ls[R\r)] (34) 

where the contribution proportional to y4[i?^(r)] is remarkably analogous to 
the Hartree energy, however, free from the self-interaction and the XCnis 
contributions are the state specific exchange-correlation functionals. They 
can be treated according to the variational principle (in some analogy with 
Ref. [37]) this is going to yield some integrodifferential equations for the 
functions R{r). This option will be considered in details elsewhere. Here we 
notice that assuming the model Slater orbital form for the functions R{r) in 
the (i-shell: 

= ^^r"-^xp(-Cr) 
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allows one to evaluate the integrals in eq. (133]) thus leading to the linear 
dependence of the latter on the orbital exponent (: 



pth ^ ^^^^ ^ . (35) 
49-76800^ ^ ^ 

^ 441-9216^ 
(here n = 3). The ratio of the theoretical values 

th/^th ^ ^093-5 91-9 ^ 

^ ' ^ 49 -76800^ 441 -9216 ^ ' 

can be compared with one extracted from the correlated and non-correlated 
experimental estimates of the parameters F2^'^ and F^^^ eqs. ( I32i) . (I3T1) which 
indicates that for some reasons the correlated model for the energies better 
agrees with the Slater model for the radial density distribution. Since as 
we mentioned the density is spherically symmetric for whatever of the states 
listed in the above Tables the only parameter characterizing the density is 
the orbital exponent (, provided the said multiplets are constructed on the 
Slater radial orbitals. In view of the linear dependency of Fj:^ on ( the state 
specific expressions for the energies and energy differences in Tables 1 and 2 
become linear functions of ( as well. 

The excitation energy expressions can be converted to the full scale den- 
sity functionals for the d-shell if one complements the above electron inter- 
action energies by the one-electron terms for the ci-shell with six electrons in 
it. The one-electron terms are (i) the kinetic energy per electron: 

(ii) the potential energy of attraction to the nucleus per electron where the 
3 in the denominator stands for the principal quantum number of the ci-shell 
under consideration: 

- f C (38) 

The average electron-electron interaction value common for all electronic 
terms is proportional to the Racah A parameter whose expression in terms 
of Fq^ and F^^ is given in the footnote to Table 1. For the iron(II) ion we 
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can set Z = 8 and the number of (i-electrons = 6 to take care about the 
core screening, then the expression for the energy becomes: 

143 

3C'-16C+15— C. (39) 
57d 

The value of ( comes then as one providing the minimum to the above func- 
tional, so that: 

in a remarkable correspondence with the Slater rules yielding for this setting 
the value of 2.08 simply by ascribing the screening increment of 0.35 to each 
electron (except one) in the (i-shell. The screening increment coming form 
the formula for A amounts 0.387. 

Including further contributions for the electron-electron interaction en- 
ergy which are now state specific yields for the ground state: 

C ^ 2.0621 (41) 

On the other hand taking one of the higher excited states whose energy 
is about 30000 cm~^ above the ground state gives: 

C ^ 2.0533 (42) 

From these estimates one can derive the following conclusion: The orbital 
exponent and thus the radial density is very weakly sensitive to whatever 
correlations. This finding is in agreement both with the accepted concept of 
correlation which attributes it exclusively to the cumulant of the two-electron 
density matrix so that there is no need to reload its manifestations on the 
density as well as with numerous demonstrations of no relation between the 
one-electron density and the correlations known in the literature (see e.g. 
Ref. [38]). 

Further analysis can be based on the observation that inserting the the- 
oretical definitions for the Fg*^ and F^^ parameters eq. (155]) into expressions 
for the excitation energies result in linear models for these energies with the 
single fitting parameter (. Two such models can be constructed: the non- 
correlated which uses only the average energies of multiple states with equal 
L and S and the correlated one which covers all ten available excitation 
energies. Fitting the excitation energies to the non-correlated model yields 
the value of 2.4823 for (. The quality of fitting with only one parameter is 
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certainly somewhat worse than that using two independent parameters F2^^ 
and F^^^ and the mean square deviation becomes 1058.0 cm~^ for the set of 
average energies of the multiplets (non-correlated fit). The value of ( which 
comes from the linear fitting procedure with the correlated energy expressions 
is 2.4604 and the mean square deviation is 1013.3 cm~^. We see that also 
in this case the correlations only marginally affect the one-electron density 
distribution and that despite some deterioration of the precision as compared 
with the two-parameter models the overall quality of the fit is surprisingly 
good. 

5.2 Summary 

Let us summarize the findings of this Section. We managed to obtain simple 
expressions for the energies of the nontrivially correlated ionic states (these 
expressions include non-dynamic correlation through the square root terms) 
with definite values of L and S as functions of a single parameter ( - the 
Slater orbital exponent for the d-shell. In the context of the accepted model 
it is the only quantity characterizing the density in the d-shell. In a sense 
there is one-to-one correspondence between the electron density of the d-shell 
and C thus the expressions for the energy can be considered as state specific 
energy functionals of the form: 

nd^-nd-C + + XCnLsindX) 43) 

In 2 5/0 

written in terms of the orbital exponent ( uniquely related to the density 
within the model used and where XCnLsi^id, C) stands for state dependent 
exchange- correlation terms as obtained by inserting the expressions for the 
Slater- Condon parameters F^^ eq. (13^ in the expressions given in Table 1 
or analogous expressions for other d-shell fillings Ref. [34]. 

As one can see our estimates of the characteristic quantity ( yield the 
values which fall into two classes depending on the type of the estimate: 
those coming from the variational estimate for the total energy of each re- 
spective state give the values close to C = 2.08 coming from the Slater rules. 
The estimates based on fitting the excitation energies to ( yield much larger 
value (much less diffuse d-shell) about 2.5 with with extremely weak influ- 
ence of electron correlation on the estimates of either of these types. These 
numerical results must be compared with other (empirical) values of the or- 
bital exponents. These, however, demonstrate a wide range of values. E.g. 
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Ref. [39] report the value C = 3.7266; Ref. [40] suggests C = 3.152; Ref. [41] 
gives C = 2.722; and Ref. [42] provide ( = 3.15 basically repeating the value 
of Ref. [40]. This indicates that either the correlated or non-correlated es- 
timates, coming from the excitation energies only, fall in the range defined 
by the Slater rules and other semi-empirical estimates. Although, the devi- 
ations between the density parameter estimates coming from different types 
of procedures also are expectable (we remind the existence of distinct ther- 
mochemical and spectral semi-empirical parameterizations) the true source 
of observed deviations is of certain interest. 

6 Range- separated treatment of electronic 
Coulomb interaction in atomic (i-shells 

Based on the idea that the short-range behavior of the e-e interactions can 
be efficiently transferred from the homogeneous e-gas to arbitrary many- 
electron systems, while the long-range e-e interactions being much more sys- 
tem specific (less transferable), Savin and Stoll suggested a generalization 
of the Kohn-Sham theory by splitting explicitly the short- and long-range 
e-e interactions [43-46]. The non-transferable long-range interactions can 
be assimilated to a wave function treatment, just like the kinetic energy in 
the conventional Kohn-Sham model, which results in a replacement of the 
non- interacting Kohn-Sham reference system by a "long-range-interacting" 
one. While in conventional KS theory the effective KS Hamiltonian has an 
exact single-determinant solution, the generalized, range-separated variant 
includes a certain amount of explicit e-e interaction and the correspond- 
ing effective Schrodingcr equation has to be solved in a multi-determinant 
form [47]. However, due to the nonsingular nature of the Ir Coulomb opera- 
tor, the solution can be converged considerably faster in both the one-electron 
and many-electron basis. 

Recent works on the range-separated hybrid methods were mostly based 
on the the separation of the Coulomb potential into short- and long-range 



30 



parts was performed according to: 




1 




(44) 



'^12/; 
1 



ri2 
erfc(//ri2) 
ri2 



erf (x) + erfc(a;) 



The treatment of the long-range exchange has been done in the Hartree- 
Fock framework, while the correlation could be treated by MP2 [48] or 
CCSD(T) [49] level, leading to a successful description of London disper- 
sion forces in vdW complexes [50] or by MCSCF level [51] to treat typical 
non-dynamic correlation problems, like the case of the H2 dissociation. A 
simpler model, where long-and short-range correlations are both handled by 
density functional approximations (RSHX - exchange-only range separated 
hybrid [52], like LC-c<jPBE of Scuseria [53]) has been recently shown to be 
quite successful in predicting magnetic coupling constants in transition metal 
systems [54]. 

In the following, we examine the behavior of the range-separated approach 
on the simple Fe(II) ion model system. 

6.1 Range separated hybrid approach 

In order to make easier the evaluation of analytical integrals and obtain 
the parameters, we decided to employ the "Yukawa" -like separation as 
proposed in Ref. [55]: 



The value of /3 — > corresponds to the absence of the long-range part. By 
contrast /3 — > 00 corresponds to the evanescence of the short-range part. The 
reach of the short-range interactions is roughly inversely proportional to the 
value of P measured in inverse bohr units. 




(45) 
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The initial assumption is that only the long-range part of the Coulomb 
interaction contributes to the non-dynamical correlations in the (i-shells so 

that only the matrix elements of f — | must be taken into account when 

the CI matrices describe the nontrivial correlation in the ci-shells. In order 
to check this assumption we have performed the following. With use of 
analytical results of Refs. [56,57] to get for the Yukawa potential the following 
expansion: 

00 I 

ki - r2r"^ exp(-/3 |ri - r2|) = 47r ^ ^ (r<r>)~2 (45) 

/=0 m=-l 

k^. (/?r<)K,^.(/?r>)F-(^) yr (^) 

where and K^^i are the modified Bessel functions of the half-integer 
index, r< = min(ri,r2), r> = max(ri,r2), and the vectors r< and r> are 
assigned correspondingly. The above expression must be inserted in the 
definition of the matrix elements of the electron-electron interaction (see e.g. 
Ref. [34]) which due to the spherical symmetry of the Yukawa potential allows 
us to express these latter in terms of the short range analogs of the Slater- 
Condon parameters. For the 3c? Slater orbitals with the orbital exponent C 
the estimates for the short range Fg'^'' and and long-range fI^ and F4 
contributions to the F^ parameters eq. ([3SD are: 

Fi^) = F2/(^);F2«=F2/«;/« + /(^) = l 
F^^ = = F4/);/) + ^;(^) = l 

= 1 — ^ — (1575/?^° + 37800/5\ + 413420/5^C^ + 

2714880/?^C^ + 11850720/3^^ + 35848960/3^C^ + 75603840/?^C^ + 
107827200/?^C^ + 94591744/3^C^ + 38578176/9C^ + 6429696C^°) (47) 

= 1 - ^^^12 (63/3'° + 1512/3\ + 16380/3^C' + 104832/3"C' + 

433888/3^^ + 1188096/3^C^ + 2101632/3^ + 2263040/3^C^ + 
1487616/?2C*^ + 559104/3C^ + 93184C^°) 
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Introducing the new variable 



t = (48) 



wc get somewhat simpler expressions for the long-range scaling coefficients 



256ii2 2304tii 8832tio 18432t9 , 

= \ \ 240^^ 

^ 91 91 91 91 

7 _ 144^5 30^=^ 15t2 



+i44t' - let'' + — + 



7 7 7 7 

(1) _ 6400^12 75520tii _ 134528t^° 1285120^^ 

^ ~ 6279 ^ 6279 2093 ^ 6279 ^ ' 

384880^^ 14160^^ 23600^^ 2451808^^ 

897 ^^3 39~ ^ 6279 

23150^^ 6910i=^ 3455^^ 
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Numerical optimization of the sum of square deviations, where the theoretical 
values are obtained under the condition that the parameters F2 and F4 under 
the square roots are respectively replaced by the long range contributions 
and ^4''*, with respect to C, and t results in the values: 

C = 2.46425 

t = 0.725436 (50) 
p = 13.0218 

The range separation parameter P is obtained by inverting the definition of 
t. These values correspond to the following scahng parameters: 

= 0.965879 (51) 
^« = 0.912733 

The precision of this estimate can be characterized as previously by the mean 
square deviation which amounts to 996 cm~^. Taking into account that the 
short-range e-e potential corresponding to (3 — 13.02 falls down to a neghgibly 
small value, say 0.001, for r = 0.5 bohr, it can be concluded that electron 
repulsion at shorter than 0.5 bohr direct space distance has an insignificant 
effect on the multiplet structure. 
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By contrast, if the theoretical values are obtained under the condition that 
the parameters F2 and F4 under the square roots are respectively replaced 
by the short range contributions F2 and , the optimization of the sum 
of square deviations with respect to ( and t results in the values: 



These values correspond to the magnitudes of the scaling parameters: 



Incidentally, the precision of the procedure singling the short-range part char- 
acterized as previously by the mean square deviation yields the value 985 
cm~^, quite similar to the long-range estimate. The reach of the "short- 
range" interactions, measured by analogous criteria as before (falling off the 
short-range Coulomb potential below 0.001) is about r = 16 bohr, which en- 
globes practically the full range for the significant densities of the d-electrons. 
It means that the long-range "tail" of the electron-electron interactions is es- 
sential to recover the correct muliplet structure. Furthermore, one can see 
that the renormalization of the F functions is in the order of 1%, confirming 
that the use of the optimal f3 implies the involvement of practically the full 
range of interactions (cf. previous Section). 

A further lesson drawn from this simple model study is that the short 
range/long range separation of the Coulomb potential is not sensitive to the 
correlations as well: the characteristic parameter of the density distribution 
C in all cases equals to 2.46 with variations in the third digit after the deci- 
mal point. Thus the short range/long range separation does not lift thus the 
strong contradiction between the estimates of the orbital exponent by the 
Slater rules or variationally from the state specific functionals eq. (jlSjl and 
those from linear fit for the excitation energies. Thunkable way out looks out 
twofold: First, the Slater rules can be thought to overestimate the screening 
(for the (i-shell the screening by the inner shells is treated to be complete, 
which yields the value 8 for the effective charge) thus leading to the values 
of ( too small as compared to those extracted from the fitting of the experi- 
mental data on excitation energies. Second, one can think that the value of 



C = 2.4653 
t = 0.0482271 
p = 0.249838 



(52) 



= 0.982441 
^(^) = 0.994545 



(53) 
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parameter Fq is for some reason much stronger renormalized as compared to 
its theoretical value eq. fl35p than those of the parameters F2 and F4. If we 
apply long-short range separation and calculate Fq at the values of (3 and C 
eq. fl52|) extracted from the fitting of excitation energies with the short-range 
parts F2 and -F4 under the square roots the fraction of the short range part 
in Fq amounts h^^^ = 0.688608 of the latter. Now if we assume that the short 
range part for some reason renormalizes to zero and thus only the long range 
part of the e-e potetial contributes to the real value of Fq then the variational 
estimate of the orbital exponent reads: 

^ = ^ _ 5 _ 143 ^ 

^ 3 2 576 ^ ' ^ ^ 

which shows some reasonable consistency with the values extracted from 
analysis of experimental spactra. Of course, this may well be a pure co- 
incidence, but possible consequences of the above hypothesis on the way 
of renormalization of the Slater- Condon parameters will be considered else- 
where. 



7 Conclusion 

In the present paper we discussed a few possible ways of avoiding the dead- 
locks of the pragmatic methods of molecular electronic structure theory based 
on the DFT, which appear due to the non- sensitivity of the basic quantity of 
the DFT - the one-electron density - to the differences in the spin (permu- 
tational) or/and spatial symmetry of the underlying many-electronic states. 
This non-sensitivity is reflected by two theorems (recent Theorems 1 and 2 of 
Ref. [14]) which formalize two basically known facts that (i) the one-electron 
density does not depend on the total spin of the many-electron state, and, 
that (ii) the one-electron density in a many-electronic state, which transforms 
according to any irreducible representation of the group acting on the spatial 
coordinates of electrons (5*0(3), 5*0(2), or their point subgroup), transforms 
according the fully symmetric irreducible representation of the corresponding 
group. 

These theorems imply that necessarily the information concerning the 
symmetry of the respective many-electronic states at hand is to be intro- 
duced into any DFT-based treatment extraneously. When it goes about the 
total spin (or equivalently about the permutational symmetry) of a many- 
electron state, we suggest to use state-specific functionals labeled by the 
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Young tableaux Tv (the rows of the irreducible representations of the uni- 
tary group U{N)) and to develop a procedure analogous to ROKS for each 
of them. In the particular case of multiple states sharing the same L and 
S in the (i-shells of transition metal ions we suggest state-specific correlated 
functionals of the density and their model based on the assumption of a 
Slater orbital form of the radial density distribution. This procedure reduces 
the functional-type density dependence to function-type dependence on the 
orbital exponent. With the use of these expressions the excitation energies 
of the many-electron states of the Fe^+ ion are reproduced with remarkable 
accuracy. The variational treatment of the proposed functionals reproduces 
with similar precision the values of the orbital exponent of the Fe^"*" ion 
prescribed by the Slater rules. Nevertheless, the estimates of the orbital 
exponent coming from the variational principle and from the fit of the ex- 
citation energies differ significantly although they fall in the range provided 
by different semi-empirical estimates. Some ideas related to conciliation of 
these two groups of estimates have been derived from analysis of the short 
range/long range separation of the electron-electron interaction potential. 
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Appendix. Permutation symmetry of the spa- 
tial function 

Since the Hamiltonian does not depend on spin variables one may wonder 
why the total spin at all affects the energy. The answer lays in the symmetry 
of many-electron wave functions with respect to permutations of coordinates 
Xi of all N electrons of the system (group Sn)- The correct wave function 
must be antisymmetric with respect to them (Pauli principle for fermions). 
This simple statement applies when the complete electronic coordinates Xi = 
(rj, Sj) ; i = 1 iV are taken as arguments of the wave function. Due to the 
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fact that the nonrelativistic electronic Hamiltonian does not depend on the 
spin projections Sj the wave function of electrons can be represented as a 
product of the spatial and spin parts dependent respectively on the spatial 
(rj) and spin (sj) coordinates only (see e.g. eq. ffTOj) ) with the antisymmetry 
requirement applicable to the entire products. In order to calculate the 
energy it is enough to know only the spatial part (multiplier) of the A^- 
electronic wave function. As one can see in eq. flTUl) the spatial parts of the 
triplet and the singlet are respectively antisymmetric and symmetric with 
respect to permutations of the spatial coordinates ri and r2 and namely this 
difference is the only real source of the differences in the energy. 

In contrast with the simple permutation symmetry properties of the com- 
plete wave functions those of the spatial multipliers are in general case some- 
what more involved. The permutation properties are conveniently described 
in terms of the Young patterns and Young tableaux. Generally the Young 
patterns are shapes formed by N boxes arranged in rows of non-increasing 
length: 



T = 



These shapes label the irreducible representations of the group S^- The 
Young pattern corresponding to the totally antisymmetric wave function 
contains only one column of the height A^. The fact that there is only one 
possibility to fill this column by electron labels from 1 to A^ corresponds to 
the one-dimensionality of the antisymmetric representation of the Sn group. 
Separation of the antisymmetric function into spatial and spin parts pre- 
defines their respective permutation properties: they must belong to the 
adjoint representations of the Sn group, since the product of two functions 
belonging to adjoint representations yields the required antisymmetric func- 
tion. The Young patterns corresponding to adjoint representations of are 
connected by 180° rotation around the bissectriss of their common upper left 
corner. Since the Young patterns which can be used for constructing the spin 
functions may contain no more than two rows, those usable for constructing 
the electronic spatial functions respectively cannot contain more than two 
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columns: 



The most remarkable feature of the Young patterns as applied to electronic 
wave functions is that they are in a one-to-one correspondence with the total 
spin, namely: the length of the one-column part of the spatial Young pattern 
equals to 25'. This allowed F.A. Matsen Ref. [23] yet many years ago to 
suggest to avoid any remark concerning the spin in (nonrelativistic) quan- 
tum chemistry context and to replace it by referencing to the permutational 
symmetry of the corresponding states. Although, it is, of course, a matter of 
terminology, within such a formulation the triplet component with Sz — +1 
would not ever arise by this hindering any possible confusion. 

The same tools can be used to describe the irreducible representations 
of the group of unitary matrices. The corresponding construct evolves as 
follows Ref. [18]: for any number of spatial orbitals M the group U{M) of 
the unitary M x M matrices acts as a "dynamical" group by transforming 
orbitals. Any given number of electrons N and any value of the total spin S 
conforming with two previous values produces an irreducible representation T 
of the group U{M). As in the case of the group irreducible representations 
of the group U (M) are labelled by the Young patterns, but the meaning of 
their elements is different. The representation by A'"-electron spatial functions 
has the tensor rank N and the corresponding Young pattern contains N boxes 
arranged in no more than two columns each of the heigth not larger than M, 
such that the first column is by 25* boxes longer than the second one. This 
irreducible representation is degenerate and its rows v can be numbered by 
distributing M orbital symbols in the above N boxes in such a way that 
they do not decrease (some ordering is assumed among them) along the rows 
and strictly increase in each column. Under this rule some orbital symbols 
in principle may appear no more than twice in a two-column pattern by this 
representing a doubly occupied spatial orbital, those appearing once represent 
singly occupied orbitals. Thus constructed Young tableaux represent states 
transforming according to the rows v of the representation T. The Young 
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tableau characterizes first of all the permutation symmetry of the state in 
that sense that the spatial part of the many electron function described by the 
Young tableau Tv is derived from the product of orbitals where each enters as 
many times as it appears in the tableau by applying the symmetrization over 
rows of the tableau and antisymmetrization over its columns. This construct 
is also known as immanant wave functions Ref. [27]. 
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Table 1: The energy expressions of the many-electron states in the d-shell of 
the Fe^^ ion. Eq is given by the expression: Eq = ridT -\-nd{nd — 1)^/2 where 
T is the kinetic energy per electron, A — Fq — 49^4 and rid is the number of 
electrons in the d-shell. 
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Table 2: The excitation energy expressions and their values for the many- 
electron states in the d-shell of the Fe^+ ion. 
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I6.F2 + I82.5F4 + 
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20300.8 
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Table 3: The average multiplct energies in the d-shell of the Fe^+ ion. 
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Table 4: The sphttings of the multiplets with coinciding L and 5" in the 
d-shell of the Fe^'^ ion. All values in cm~^. 
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